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Resonating valence bond (RVB) states are of crucial importance in our intuitive understanding 
of quantum spin liquids in 2D. We systematically classify short-range bosonic RVB states into sym- 
metric or nematic spin liquids by examining their flux patterns. We further map short-range bosonic 
RVB states into projected BCS wave functions, on which we perform large-scale Monte Carlo simula- 
tions without the minus sign problem. Our results clearly show that both spin and dimer correlations 
decay exponentially in all the short-range frustrated (non-bipartite or Z2) bosonic RVB states we 
studied, indicating that they are gapped Z2 quantum spin liquids. Generically, we conjecture that 
all short-range frustrated bosonic RVB states in 2D have only short-range correlations. 



Introduction: Quantum spin liquids are exotic in- 
sulators which cannot be adiabatically connected into 
a band insulator and which can support fractional- 
ized excitations pQ. Introduced by Anderson nearly four 
decades ago [2 , the resonating valence bond (RVB) state 
on the triangular lattice is the first example of quantum 
spin liquids in more than one dimension. Since then, 
there has been keen interest in searching for such exotic 
states of matter in real materials as well as in micro- 
scopic models, especially after exciting connections be- 
tween quantum spin liquids and the mechanism of high 
temperature superconductivity were suggested [3]-[5]. 

Recently, there has been a surge of numerical simu- 
lations on simple models reporting convincing evidence 
of the existence of fully gapped spin liquids [BUTT], all of 
which are believed to be in the same class of short-range 
bosonic RVB states. Nonetheless, the nature of short- 
range bosonic RVB states on various frustrated lattices 
has not been explicitly revealed [T2] [T3], mainly because 
of the so-called minus sign problem in Monte Carlo (MC) 
simulations of those bosonic short-range RVB states with 
frustration. Because of their conceptual importance in 
pictorially understanding quantum spin liquids and their 
direct relevance in recent numerical simulations, it is 
highly desired to unambiguously demonstrate the nature 
of these short-range bosonic RVB states. 

In this Letter, we systematically classify short-range 
bosonic RVB states by examining their flux patterns [14] . 
For instance, for the Kagome lattice we establish that 
there are only four symmetric RVB states when consid- 
ering only nearest-neighbor (NN) valence bonds, as is 
shown in Fig. 1. Then, we show that these bosonic short- 
range RVB states (including a class of RVB states with 
valence bonds longer than NN) can be exactly mapped 
into projected BCS wave functions [T5HT7] on which we 
perform large-scale MC simulations without the minus 
sign. For frustrated short-range RVB states, our simula- 
tions on corresponding projected BCS states convincingly 
show that both their spin and dimer correlations decay 



exponentially, indicating that they are fully gapped Z2 
spin liquids pi IT9]. 

Bosonic RVB states: We consider the following 
bosonic RVB states with NN and possibly next nearest- 
neighbor (NNN) valence bonds 

h/>RVB> = Ei c >> \c) = (-v Se II Aily). (!) 

c (ij)ec 

where c labels valence bond configurations and 8 C repre- 
sents the number of bond-crossings in c [the factor (— l) 6c 
is nontrivial only for RVB states with valence bonds 
beyond nearest-neighbor sites.] Here \ij) = (|tiij) — 
\^i\j))/y/2 is the spin-singlet wave function (or valence 
bond) on (ij) and we assume \fij | to respect all the lattice 
symmetries. Note that Eq. Q represents a "bosonic" 
RVB state in the sense that |tiij) = lijti)- Since 
\ij) = — \ it is sufficient to consider — —fji. The 
wave function in Eq. ([TJ possesses a gauge symmetry: 
I^rvb) is invariant, up to a phase, under the transfor- 
mation fij — )> e loti fije zaj . In the following, we focus on 
time reversal invariant RVB states for which all fij can 
be chosen as real; signs of can be represented by ori- 
ented arrows on graphs: an arrow pointing from i to j 
means that > 0, as shown in Fig. [I] We further define 
flux (j) p = 0, 7r (mod 2tt) for plaquette p through 

cc 

n sgn(/ j/c ) = exp(^ p ), (2) 

(jk)£p 

where cc means that the counterclockwise order of (jk) 
is taken in the product above and sgn is the sign func- 
tion. It is clear that (f) p is gauge invariant for even-length 
plaquette p. However, the gauge transformation with 
exp(zaj) = i on every site j changes <j) p to (j) p + ir for 
all odd-length plaquette p. In other words, the two wave 
functions with flux patterns {(j) p } and {<p p + (— l) Up 7r} 
(n p is the length of plaquette p) actually represent the 
same state [20] . 
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FIG. 1: (a) The flux patterns {4> P } of the only four symmet- 
ric NN-RVB states on the Kagome lattice. Here e x and e y 
represent the unit vectors, (b) The flux patterns in the 

corresponding projected BCS states. 



There are two questions to be answered concerning the 
wave function in Eq. ([!]). First, is it a symmetric spin liq- 
uid respecting all the symmetries of the lattice in ques- 
tion? Second, do various correlations decay in power law 
or exponentially? The first question can be answered 
by examining its flux pattern {</> p }. If the flux pattern 
{4>p} is invariant up to the addition of {( — l) Up 7r}, under 
all lattice symmetry transformations such as translations, 
rotations, and reflections, the corresponding RVB state is 
then a symmetric spin-liquid state . We label RVB states 
whose longest valence bonds are between NN (NNN) sites 
as NN-RVB (NNN-RVB) states. On the Kagome lattice, 
we identify four NN-RVB states as symmetric spin liq- 
uids, as shown in Fig. [I] On the triangular lattice, only 
two symmetric NN-RVB states are found, as shown in 
Fig. [2j On the square lattice, there are two symmet- 
ric NN-RVB spin liquids and four symmetric NNN-RVB 
states, as is shown in Fig. [3J 

For these symmetric RVB spin liquids, it is not known 
a priori whether various correlations decay in power law 
or exponentially. Generically, numerical MC simulations 
are capable of revealing those features [2 1. The correla- 
tions of a physical quantity O are given by 



(OiO s ) = 



{lpR\B\OjOj I^RVB) 

(^rvbI^rvb) 
£ cc ,(c|c') 



<c|c') 



E c ,c' < c l c '> 



(3) 



where (c\c f ) (| (c\c f ) |) can be taken as statistical weight 
in MC simulations when they are positive (negative). For 
instance, for the square lattice NN RVB state with {<p p = 
0}, (c\c f ) > for any c and c', on which large-scale loop 
algorithm MC simulations [22 were performed recently, 
reporting evidence that this RVB state is critical with 
power-law decaying dimer correlations [23 , [24] . 



Kagome NN-RVB state 


A 
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D 




0.6 


0.6 


0.6 


0.7 




1.2 


1.0 


1.0 


0.9 


E/J 


-0.393 


-0.36 


-0.357 


-0.386 



TABLE I: The spin (£ s ) and dimer correlation lengths 
of the four symmetric states on the Kagome lattice shown 
in Fig. ^ Here E labels the variational energy per site of 
those symmetric states for the Kagome NN antiferromagnetic 
Heisenberg model H — J ^2 Si • Sj. 

(ij) 



However, it is impossible to choose (c\c f ) > for all c 
and d for NN RVB states on non-bipartite lattices (e.g. 
the triangular lattice) or NNN RVB states on bipartite 
lattices (e.g. the square lattice). Such states are exam- 
ples of frustrated RVB wave functions defined as ones 
whose valence bonds form non-bipartite graphs. It is 
clear that loop-algorithm MC simulations on frustrated 
RVB states suffer from the minus sign problem in the 
variational level. In the following, we shall show that 
the RVB states in Eq. ([I]) can be exactly mapped into 
Gutz wilier projected BCS states, which are friendly to 
MC simulations without the minus sign problem. 

Projected BCS states: It was known that the vari- 
ational Monte Carlo method has been quite successful in 
simulating Gutzwiller projected BCS wave functions. We 
consider the following projected BCS wave functions: 



I^p-bcs) = Vcexp 



to) 



10} , (4) 



where c\ a are electron creation operators, |0) is the vac- 
uum, Vg is the Gutzwiller projection onto singly oc- 
cupied states, and g^ = g^ which are assumed to be 
real. A similar gauge symmetry exists for the pro- 
jected BCS wave functions: the wave function is invari- 
ant, up to a phase, under the transformations: gjk — » 
ex.p(iaj)gjk exp(za£). For time reversal invariant states 
with real g^, we define fermionic fluxes $ through 
U(jk)ep s ^ n (9jk) = exp(i^) with <^ = 0,tt (mod 2tt). 
Similarly, the flux pattern and + ( — l) np 7r} are 

equivalent through the gauge transformation exp(zaj) = 
i on all sites j. As shown in details in the Appendix, we 
obtain I^p-bcs) = |^rvb) when the following conditions 
are satisfied: 



\9jk\ = \fjk\ 



bp (j) f p 



(5) 



where p labels all possible elementary plaquettes. 

Ground state correlations: We have investigated 
a number of short-range frustrated RVB states by per- 
forming MC simulations on their corresponding projected 
BCS wave functions. In our numerical calculations, we 
mainly focus on spin as well as dimer correlations and 



(A) (B) 

FIG. 2: The flux patterns {cj) p } of the two symmetric NN-RVB 
states on the triangular lattice. 

study whether they fall off in power law or exponentially 
at large distance. The spin and dimer correlations are 
defined as follow: 

S(l) = (s fi • s rl+r ) , 

where e a labels lattice vectors. 

We have studied all the four symmetric NNN-RVB 
states on the square lattice, the four symmetric NN-RVB 
states on the Kagome lattice, and the two symmetric 
NN-RVB states on the triangular lattice. As discussed 
in details later, for all these frustrated RVB states our 
MC simulations convincingly show that their spin and 
dimer correlations fall off exponentially with correlation 
length in the order of one lattice constant, indicating 
that they are all gapped Z^ quantum spin liquids. (Note 
that different symmetric spin liquid states on the same 
lattice may be distinguished by numerically computing 
local correlations.) We conjecture that our results are 
generic: all frustrated short-range RVB states in 2D are 
fully gapped. 

The Kagome lattice: According to the flux pattern 
(either {(j)p} or {</>p}) on the Kagome lattice, it is 
straightforward to show that there are four symmet- 
ric NN-RVB spin liquids, as shown in Fig. [T] We 
have computed the spin and dimer correlations for all 
the four symmetric NN-RVB spin liquids. In Fig. |4j 
we plot the spin and dimer correlations as a func- 
tion of distance (I) in one of those NN-RVB symmetric 
states, i.e. the counterclockwise-counterclockwise NN- 
RVB state (A) shown in Fig. [I] The MC calculations 
are carried out on a lattice with 18 x 18 x 3 sites and 
with periodic boundary conditions. It is remarkable that 
the correlations decay extremely fast. From Fig. [4j it 
is clear that both spin and dimer correlations decay ex- 
ponentially with distance. The spin correlation length 
£ s for this NN-RVB state is about 0.6 lattice constants. 
The dimer correlation lengths £d,x an d id,y f° r D x (le x ) 
and D x (le y ) are nearly equal, which are about 1.2 lat- 
tice constants. The spin and dimer correlation lengths 
[£s and £ d = (£ diX + fd >3/ )/2, respectively] for all the four 
symmetric NN-RVB spin liquids are listed in Table [TJ The 
correlation lengths are all in the order of one lattice con- 
stant, indicating that they are fully gapped Z2 quantum 
spin liquids. This is consistent with the recent density 
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FIG. 3: The flux patterns {</> P } of (a) the two symmetric NN- 
RVB states and (b) the four symmetric NNN-RVB states (b) 
on the square lattice. 



square NNN-RVB state 
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D 




1.1 


0.8 


0.7 
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£,d,NN 


1.2 


1.6 


1.3 


1.4 


£cZ,NNN 


0.6 


0.6 


0.6 


0.5 


E/Jt 


-0.344 


-0.219 


-0.237 


-0.239 



TABLE II: The correlation lengths and variational energies 
of the four symmetric NNN-RVB states on the square lattice. 



matrix renormalization group evidence that the ground 
state of the Kagome Heisenberg antiferromagnet is a fully 
gapped quantum spin liquid [HJ [9] with correlation lengths 
of about one lattice spacing [9]. 

To get a sense of which of the above states is the 
best variational wave function for the Kagome antifer- 
romagnet described by H = J^^Si * Sj, we com- 
pute their variational energy per site as is shown in Ta- 
ble |TJ Among the four symmetric NN-RVB spin-liquid 
states, the counterclockwise-counterclockwise NN-RVB 
state has the lowest energy for the Kagome antiferromag- 
net, which is — 0.393J per site. This energy still differs 
from the density matrix renormalization group result, in- 
dicating that longer-range valence bonds are important 
in describing the Kagome antiferromagnet. 

The triangular lattice: It turns out that there are only 
two symmetric NN-RVB spin liquid states on the triangu- 
lar lattice, whose flux patterns {4>p} are shown in Fig. [2J 
For both states, the spin and dimer correlations decay ex- 
ponentially with distance, with £ s = 0.7 and ^ = 1.0 for 
the state (A) shown in Fig. [2] and £ s = 1.0 and ^ = 1.6 
for the state (B) shown in Fig. [2] Both NN-RVB states 
are then gapped Z2 spin liquids. The correlation lengths 
on the triangular lattice are somewhat longer than those 
on the Kagome lattice, which is expected due to more 
geometric frustrations in the Kagome lattice. 

The square lattice: It was shown recently that the un- 
frustrated NN-RVB state on the square lattice is a criti- 
cal state with power-law dimer correlations even though 
its spin excitations are gapped[23-25 . To have a fully 
gapped spin-liquid phase on this lattice, it is necessary 
to include frustration in short-range RVB states. In this 



FIG. 4: The spin and dimer correlations as a function of 
distance in the NN-RVB state (A) shown in Fig. [I] 



Letter, we consider to include NNN valence bonds, which 
is partly motivated by a recent study establishing that 
fully gapped spin-liquid ground states are realized in the 
generalized [26 quantum dimer models [27 with NN and 
NNN dimers on the square lattice. (Note that the nature 
of NNN RVB states in Eq. ^ without the factor (-l) 6c 
remains unknown due to the lack of mapping between 
them and projected BCS states.) 

From the flux pattern {</> p }, we have identified four 
symmetric NNN-RVB states on the square lattice, as 
shown in Fig. (3^b). For each of these four states, there 
is an additional parameter labeling the wave function, 
namely the ratio 7 = |/nnn//nn|- We take 7 = 1 in 
our MC simulations. For 7 = 1 both spin and dimer 
correlations decay exponentially with distance. The cor- 
relation lengths are listed in Table [TTj where ^,nn and 
£d,NNN mean the correlation lengths of NN and NNN 
dimers, respectively. At 7 = 1 (more generally 7 > 0), 
we conclude that the four symmetric NNN-RVB states 
are fully gapped Z2 spin-liquids, which is consistent with 
the recent numerical evidence of fully gapped spin-liquids 
in the J1-J2 square Heisenberg antiferromagnet [TO ] [TT] . 
The variational energies of these NNN-RVB states in 
unit of J 1 for the J1-J2 square Heisenberg model with 
J 2 — Ji/2 are shown in Table |TT| 

Nematic RVB spin-liquids: We have studied fully 
symmetric short-range RVB spin-liquids on various lat- 
tices. An interesting question is whether short-range 
RVB states could be nematic spin-liquids which are trans- 
lationally invariant but break lattice point group symme- 
try. The answer is yes. On the Kagome lattice, we iden- 
tified that there are only four NN-RVB states which are 
nematic spin-liquids, as shown in Fig. [5|a). On the tri- 
angular lattice, there are two nematic NN-RVB states, as 
shown in Fig. [5jb) . Our MC simulations show that the 
correlation functions of spins and dimers in these states 
decay exponentially with distance but the Cq v symmetry 
of both lattices is broken. They are fully gapped nematic 
spin-liquids, in contrast with gapless nematic spin-liquids 
on the triangular lattice studied in Ref. [28] . 

On the square lattice, we found six nematic NNN-RVB 
spin-liquid states, which are shown in Fig. [5|c). These 



(c) 




FIG. 5: The flux patterns {4>p} of (a) the four nematic NN- 
RVB states on the Kagome lattice, (b) the two nematic NN- 
RVB states on the triangular lattice, and (c) the six nematic 
NNN-RVB states on square lattice. 

spin-liquid states keep the translational symmetry but 
breaks the C± v rotational symmetry of the square lat- 
tice. Again, our MC simulations show that they are fully 
gapped nematic spin-liquids. 

Concluding discussions: On the square (or honey- 
comb) lattice, there are two kinds of symmetric NN-RVB 
spin-liquids. The {(j) p = 0} NN-RVB state is unfrustrated 
with power-law decaying dimer correlations [23] [24]. For 
the {(j)p = 7r} NN-RVB state, our MC simulations on the 
square lattice with 40 x 40 sites implies that its dimer 
correlation decays in power law even though we need 
further studies on systems with larger sizes to deter- 
mine the power exponent accurately. [Here the power-law 
decaying dimer correlations are expected since bipartite 
RVB states are effectively described by an emergent U(l) 
gauge theory.]. 

A recent loop-algorithm MC study shows that the un- 
frustrated NN-RVB states on the cubic and diamond 
lattices show magnetic long-range order [29]. Properties 
of NN-RVB states on 3D frustrated lattices remain un- 
known partly due to the minus sign problem in loop- 
algorithm MC simulations. We can generalize the map- 
ping between short-range bosonic RVB states and pro- 
jected BCS states discussed in the present Letter to three 
dimensions. Consequently, we can solve the minus sign 
problem for a class of frustrated short-range RVB states, 
which is a significant step towards understanding the na- 
ture of short-range frustrated RVB states in 3D. 
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Note added: After the completion of the present Letter, 
we notice Ref. [30] where similar results were obtained 
for a subset of NN RVB states classified in the present 
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Letter. 
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Supplemental Material: proof of exact mapping between bosonic RVB and projected BCS states 



Now, we shall prove the mapping between a class of short-range RVB states and Gutzwiller projected BCS wave 
functions [16]. The former is given by 



c (ij)ec 



= D-D 

c 

= ^(-l)*.(-l)Pc 

c 

= £(-l)'.(-l)Pe 



n fij( s i- s n 

(ij)ec 



Iff • • • t> 



n h 3 {s- - s-)4^ 

(ij)ec 
(ij)ec 



|o), 

10} , 



(6) 
(7) 
(8) 
(9) 



where p c is the signature of c (some ordering in (ij) for each c is implicitly assumed.) 
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Let's first consider NN-RVB states, for which S c = for any c. The NN-RVB states are given by 



NN-RVB 



(-i) Pc n u4Ai-4A t ) 

(ij)ec 



|o) = £|c>. 



(10) 



The corresponding projected BCS wave function Eq. Q with \gij\ = \ fij\ can be expanded as follows: 



I^P-BCS) = V G 



(ij) 



10), 



(ij)ec 



|o> = 5» F . 



(11) 



Notice that the second step above is a consequence of the no-double occupance projection, which excludes the possi- 
bility that one site takes part in two valence-bonds. 

Suppose | Co) = \cq)f for some Co by choosing appropriate signs of gij for (ij) G Co, which can always be done. We 
have 



(-l)^o ]J sgn(/ ii )= J] sgnfoy)- 
It is clear that two wave functions I^nn-rvb) and |^p-bcs) are equal if 

(-1)* n s s n (/^) = n s ^j) 

(ij)ec (ij)Gc 



for any c. Eq. ( 13 ) is equivalent to 



(_l)Pc+Pc "Q sg n(/ ii )= J] sgn(^), 

(ij)Gc+c (ij)^ c + c o 



(12) 



(13) 



(14) 



where c + cq represents a transition graph obtained by superposing two valence bonds configurations c and on the 



same lattice. The transition graph consists of even- length loops of valence bonds. It is clear that Eq. (14) is satisfied 
if 



n s Mdij) = - n s s n (/^) 

(ij)ep (ij)ep 



(15) 



for every even- length plaquette p. The minus sign in Eq. (15) comes from the signature (— 1)p<=+p<!o when the ordering 
of (ij) in c + Co is taken in a counterclockwise fashion. For Eq. (15) to be true, it is sufficient to have 



(16) 



for all possible elementary plaquettes p, which could be even or odd-length. Note that open boundary conditions of 
the lattice in question is implicitly assumed here. For lattices with periodic boundary conditions in both directions, 
four projected BCS wave functions are needed to form a linear superposition in order to represent a short-range RVB 
state [171 US]- Nonetheless, the boundary conditions of lattices do not qualitatively affect the large-distance correlation 
functions of local physical quantities. 

Now, we consider the case where both NN and NNN valence bonds are allowed in valence bond configurations 
and the RVB states are called NNN- RVB states. It is clear that each cross generates a twist from the configuration 
without that cross, which results in the factor (—1) from the signature ( — l) Pc . The factor (—1) from the signature 
(— l) Pc can be compensated by the factor (— l) Uc . This proves that the NNN-RVB states can also be mapped exactly 
into projected BCS wave functions with pairing between NN and NNN sites. 



